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SUMMARY 


A preliminary calculation of an airplane spin-recovery motion has 
heen made and the calculation technique has heen evaluated as a possihle 
means for Investigating spin and recovery motions for extreme fligjht 
conditions difficult or in^jossihle to simulate by free -spinning-tunnel 
tests. The recovery was calculated step by step by using modified wind- 
tunnel rotary-balance measurements^ applicable eqmtione of motion^ and 
spin-geometry relationships. Diffictilties encountered in using the 
rotary-balance data in the calculations are discussed, and it is pointed 
out that certain inconsistencies must be cleared up before the method 
can be accepted as adequate to give detailed spin-recovery motions for a 
specific airplane. Time-history plots are presented which represent 
variations, after rudder reversal, for attempted recovery, of ensuing 
attitudes, velocities, and accelerations of the airplane with respect to 
the airplane body axes, the earth, the relative wind, and the resultant 
axis of the spinning motion. The time -history ciarves indicate initial 
small oscillations of the factors from their original valTjes, with these 
oscillations increasing gradually in amplitude and with nearly all the 
large significant changes occurring near the end of the time required for 
recovery. 


INTRODUCTION 


Many investigations have been made in the Langley 20-foot free- 
spinning tunnel in which the natures of airplane spins and recoveries 
have been determined experimental ly with dynamically scaled-down models. 
The results of some of these investigations are presented in references 1 
to ' 5 . The trends in airplane design, however, have increased the diffi- 
culty of maMng dynamic -model spin tests and interpreting the results, 
either because of the nature of the motion obtained or because of tunnel 
limitations. Detenninatlon of spin and recovery motions for some current 
extreme airplane fll^t conditions cannot be made by the dynamic-model 
techrdLque as, for exanple, motions of airplanes at extr e mely hi^ alti- 
tudes at which the model wing loadings requiMd to simulate the airpla n es 
dynamically are inpracticably hl^ for testing in the free -spinning 
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tunnel. Also, the free -spinning-tunnel technique cannot he used to 
Investlga-te spinning motions which may he encountered in flight at hl^ 
Ifeich nunibers. Thus, it appears desirable to he ahle to determine the 
motion of an airplane during the spin and recovery analytically. 

The results of some previous analytical considerations of spin 
recovery were reported in reference 6 in which step-hy-step calculations 
were made to study the relative efficiencies of the three airplane con- 
trols (rudder, ailerons, and elevator) in achieving recovery from an 
estahllshed spin; however, because of a lack of more comiplete aerodynamic 
data, this work was handicapped by the necessity of using estimated aero- 
dynamic forces and moments based on the assun^tlon of constant derivatives. 
The present investigation is a preliminary evaluation of a calculation 
technique which makes vise of aerodynamic data obtained on a wind-tunnel 
rotary balance and in which the assumption of constant aerodynamic 
derivatives is not necessary. In this investigation, applicable equations 
of motion and spin-geometry relationships were used with the rotary- 
balance aerodynamic data to calculate step by step the details of a spin- 
recovery motion for an unswept-wing filter -airplane configuration. 

Because of certain inconsistencies in the rotary-balance data (similar to 
those previously discussed in ref. 7)> some modifications were made to 
the aerodynamic force and moment coefficients used in the calculations, 

Hie inconsistencies are believed to be primarily due to difficulties 
involved in mounting the airplane model in order to obtain correctly the 
aerodynamic data to simulate a steady spin. These inconsistencies and 
the modifications are discussed in detail subsequently. 

The calculated recovery motion is presented as time-history plots of 
the variations, after rudder reversal for attempted recovery, of ensvdng 
attitudes, velocities, and accelerations of the airplane with respect to 
the airplane body axes, the earth, the relative wind, and the resviltant 
axis of the spinning motion. 


SYMBOLS 


The calculated recovery motion is presented with respect to the air- 
plane body system of axes and to space and wind attitude angles (fig. l) 
and also with respect to the axis of resultant rotation of the motion. 

Cv longitudinal-force coefficient, — — 

|pVr2s 


later aJ--f or ce coefficient, — ± 

1 p 


Cy 
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m. 


•'ll 


nortnal. -force coefficient, -j- 
rolling-moment coefficient. 


JpVr2s 
L 


JpVR^Sb 


pitcMng-moment coefficient. 


M 




yawing-moment coefficient. 


N 


IpVE^Sb 


X longitudinal force, lb 

Y lateral force, lb 

Z normal force, lb 

F resiiltant force, lb 

L rolling moment, ft-lb 

M pitching moment, ft-lb 

N yawing moment, ft-lh 

S wing area, sq, ft 

b wing span, ft 

p air density, slugs/cu ft 

u,v,w components of resiiltant linear velocity Vg along the X, Y, 

and Z body axes, respectively, positive in directions of 
positive X-, Y-, and Z-forces, ft/sec 

Vj^ resultant linear velocity of airplane, ft/sec 

p,q,r components of resultant angular velocity ft about the X, Y, 

and Z body axes, respectively, positive in directions of 
positive L-, M-, and N-moments, radians/sec 

ft resultant angu 1 ar ^ velocity (if axis of resultant rotation is 

vertical, ft = radlans/sec 
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IJ- 

m 


c 


c 


x/5 


z/c 




airplane relative-densily coefficient, 

mass of airplane, slugs 

S 


m 

pSb 


local chord, ft 

mean aerodynamic chord, ft 

ratio of distance of center of gravity rearward of leading 
edge of mean aerodynamic chord to mean aerodynamic chord 

ratio of distance between center of gravity and X body axis 
to mean aerodynamic chord (positive when center of gravity 
is below X body axis) 

radii of gyration about X, Y, and Z body axes, respectively, ft 


moments of inertia about X, Y, and Z body axes, respectively, 
SlTlg-ft^ 



% - iz 

nOb^ 


inertia yawing-moment parameter 


inertia rolling-moment parameter 


^Z - % 

mb^ 

S 


^e 


inertia pitching-moment parameter 


acceleration due to gravUy, talsen as 52 . 1 ? ft/sec^ 

vertical component of total angular deflection of X body axis 
from a reference position in horizontal plane, positive 
when sense of deflection is airplane nose Toward from 
reference position, deg or radians as indicated (fig. 1(a)) 

horizontal ccmponent of total angular deflection of X body 
axis from reference position in horizontal plane, positive 
when clockwise as viewed from vertically above airpla n e, 
deg or radians as Indicated (fig. 1(b)) 
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ai 


“5 

“i 


a 


P 



8Vr 


a 


component in YZ "body plane of angular deflection of Y "body 
axis from horizontal plane, positive when clockwise as 
viewed from rear of airplane (if X body axis is vertical, 
0 g is measured in horizontal plane) , deg or radians as 

indicated (fig. l(c)) 

angle between Y body axis and horizontal, positive for erect 
spins when right wing is down and for inverted spins when 
left wing is down, deg 


acute angle between Y body axis and a plane perpendicular to 
axis of resiiltant rotation, sin ^ deg 


acute angle between X body axis and vertical, deg 


acute angle between Z body axis and horizontal, deg 

angle between X body stxis and axis of resultant rotation, 

cos“^ — or sin“^ deg 

ft ft ' 

angle of attack between relative wind Vp projected into 
XZ plane of symmetry and X body axis, positive when relative 
wind comes dfrom below XY body plane, deg 

angle of sideslip at center of gravity (angle between relative 
wind Vj and XZ plane of symmetry), positive when relative 
wind comes from ri^t side of plane of symmetry, deg 

angle of Inclination of relative wind Vr from XY body plane, 

positive when inclined below XY body plane, sin"^ ~, deg 

Yr 

angle between relative wind Vr projected into XY body 
plane and X body axis, positive when Vr comes from 

right side of plane of synmetry, sin”^ deg 

Yr 

angle of inclination of axis of resultant rotation from 
XY body plane, positive when inclined below XY body plane, 

sin”^ deg 
ft^ 
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5a 


Pi 

t 

At 

R 


a 


V 


•'a 

^Vr 


angle "between axis of restiltant rotation projected into 
XT "body plane and X tody axis, positive when vector of 
axis of resultant rotation lies on right side of 

plane of symnetry, sin"^ -====, deg 

\|PT7 

angle of Inclination of resultant force vector from XT hody 
plane, positive when inclined helow XX hody plane, deg 

angle between restoltant force vector projected into XI body 
plane and X body axis, positive when resultant force 
vector comes from right side of plane of symmetry, deg 

time, sec 

time Increment, sec 

spin radius in developed spin, distance from spin axis to 
center of gravity, ft 

angle between direction of resultant linear velocity and 
axis of resultant rotation (helix angle in developed spin), 
deg 

velocity of airplane center of gravity along axis of resultant 
rotation, Vr cos a, ft/sec 

angle between axis of resultant rotation and vertical, deg 
angle between resviltant linear velocity and vertical, deg 


I,U,III,IV,V auxiliary angles in spherical triangles used in determining 
a (fig. 6) and kq (appendix B), deg 

Vy vertical conponent of velocity of airplane center of 

gravity, ft/sec 

Subscripts: 

n indicates end of time Increment 

n-1 indicates beginning of time increment 

A dot over a symbol represents derivatives with respect to time; for 

example, d = — . 

dt 
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The attitude angles 0 q, \lre> and (Euler's space angles) as 
defined iised herein are Intended to provide a means of interpreting 
any changes from the reference attitudes for either an initially erect 
or an initially inverted spin. Under this general system, the reference 
attitude for is arhitrarily taJcen as zero for either an erect or an 

inverted spinj the airplane would go from an erect to an inverted atti- 
tude or vice versa if Gg hecomes greater negatively than -90° or if 
exceeds +90°. 


GENERAL CQNSIDERATIQNS 


A spin is, in gene2ral, a motion in which an airplane in flight 
descends rapidly toward the earth while rotating with angijlar velocity 
ijrg about a vertical axis and while at an angle of attack above that at 
which the stall occurs . The spin, depending upon the attitude of the 
airplane in space, may be erect or inverted. If the spin is erect, the 
angle of attack is positive and, if the spin is inverted, the angle of 
attack is negative. The airplane must be designed so that recovery from 
a spin (either erect or Inverted) can be achieved following prescribed 
manipulations of the airplane controls to enable the pilot to regain 
normal control of the airplane. OJhe motion of a specific airplane during 
a recovery is apparently dependent on the spin motion before the controls 
are moved, on the particular control manipulation en^iloyed for recovery, 
and on the effectiveness of these controls (refs . 2 to 5) • 

Based on dynamic tests on many models in the Langley free-spinning 
timnel, it appears that the motion following control movement for 
recovery may be extremely varied. The variations may range from a motion 
in which the model, immediately after control movement, decreases its 
rate of rotation and noses down rapidly for a recovery to one In which 
movement of the controls has little or no effect on the spinning motion. 
In some instances, after controls are moved for recovery, oscillations 
may occur in the attitvide of the model and in its rate of rotation, with 
an over-aJ-1 trend toward a decrease in both rate of rotation and angle of 
attack which eventiially evolves into a steep low-angle-of -attack pull-out 
dive without rotation (ijre becomes zero) . Recovery, however, is sometimes 
considered as having been achieved if the model appears to have entered a 
spiral glide or an aileron roll at an angle of attack below the stall. 

In some cases, recovery from the spin may be associated with a quick 
rolling or pitching motion to an inverted attitxide, and recovery is con- 
sidered to have been achieved from the original erect spin regardless of 
whether \Irg has changed its direction in space. From a practical stand- 
point it would be necessary that the pilot of a corresponding airplane be 
able to exercise sirfficient care to avoid entering another spin after 
having achieved recovery from an original spin. These general considera- 
tions would also apply were the original spin inverted rather than erect. 
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Wlien a free-spinnirig-’turinel model has recovered from a spin within 
2^ turns after control manipulation, the airplane being simulated hy the 

model is customarily considered to have satisfactory spin-recovery char- 
acteristics. This figure has been deduced from escperience and is based 
on comparisons between free-spinning-tunnel test results and available 
full -scale-airplane spin data. 


METHODS AMD CALCULATIONS 


The airplane is considered as a rigid body in space having six degrees 
of freedom. Equations of motion are used which relate the linear and 
au gul ar velocities of the airplane along and about the airplane body axes 
to the attitvide of the airplane with respect to the earth. The equations 
are equivalent to equations given in reference 8 and in various dynamics 
texts, except that no product of Inertia terms appear because the airplane 
body axes are assumed to be the principal axes. The equations are 


u = Gy - g sin 0_ + vr - wq 

2pb ® 


V 

^ ” 2pb 


Cy + g cos 9 g sin +. wp - lu: 


Vr2 

❖ = Cg + g cos Qg cos Pg + uq. - vp 


i Cl + 


% ~ 


qr 


‘•X 


4 = -^Cm + -'^ 


± = 


2|ikY 


_ ’'E 


C„ + 


Ix-% 


pq. 


2pkzi 


The airplane is considered to be initially in an erect steady spin 
turning to the pilot's ri^t in space (positive 1^^) . Since in a steady motion 
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all time derivatives of the variables are zero, the equations of motion 
show that 0Q and must he constant dirring a steady spin. Therefore, 

inasmuch as 0g and are zero, all the spin rotation in terms of 

Euler's angles must he due to ij’g, which is a rotation about a vertical 

axis. The conponents of the total angular velocity about airplane body 
axes are given in terms of the rates of change of the Euler angles by 
the relationships (ref. 9) 

P = ^e - ^e ®e 
' q = 0Q cos + tg cos 0g sin 

r = -0g sin 0g + tg cos 0g cos 

These relationships for p, q, and r and the following relationships 
were used to find values for and for the velocity factors in the 

equations of motion when time equaled zero: 

0g = sin“^ - (see fig. 3) 

® cos 0e 

Vr =^v/ + (Riirg)^ 

V = Vjj sin p 


u = Vj^ cos p cos a 


and 


w = u tan a = Vj^ cos 0 sin a 


where 


R « 


g tan 0^ 

^e 


p w 0 - cr 
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a, « 90 - 0 g 


a *» taaa 


-1 Bte 


and where 0 g, 0 , Vy, and ijrg were obtained from available free- 

splnning-tunnel test results (avereige values) of a l/20-scale dynamic 
model of the airplane being considered. 

l&iss characteristics and control settings for the airplane and some 
of the Initial erect-spin characteristics are presented In table I. Tha 
airplane control manipulation simulated for recovery In the calculations 
Is reversal of the rvidder from full with (pro-spin) to full, against the 
spin, with no further control adjustments being simulated. 

IMpubllshed aerodynamic data were available for a l/lO-scale model 
of the configuration being considered from rotary-balance measurements 
made in the Langley 20 -foot free-spinning tunnel. A drawing of this model 
is shown in figure 2 . The rotary balance used consists of a six-component 
strain gage that measures forces along the three airplane-model body axes 
and moments about these three axes while the model Is being forcibly 
rotated at a constant rate about a vertical axis in the vertical airstream 
of the timnel. The angles of attack and sideslip of the rotary-balance 
model may he set at angles from 0 ° to 56 (P. A con5)lete description of 
this rotary balance is contained in reference 7- As obtained, the avail- 
able rotary-balance data incl\ided values of the force and moment coeffi- 
cients for wide ranges of p, a, and fl with the model rudder neutral 
and also values of the force and moment coefficients for rudder-with and 
rudder-against conditions at cP sideslip for the ranges of a and ft. 

The control settings used on the model were somewhat different from those 
used in the free-spinning dynamic-model tests mentioned previously 
(table 1) . The model elevator was set at 50° up and the ailerons were 
set at 120° against the spin (in corresponding airplane, stick left in 
spin to pilot’s right). All these data were obtained from rotation tests 
in which the model was rotated about a vertical axis through its center 
of gravity (zero spin radius) . 

The aerodynamic data obtained include no effects of accelerations 
such as are accounted for by the dot stability derivatives in studies of 
airplane stability or motion. The effect of the lack of these accelera- 
tions in obtaining the aerodynamic data is probably negligible when 
attenpting to simulate a fairly steady spin, and although the effects on 
a calculated spin-recovery motion are not known, they may not be appreci- 
able. In obtaining data with the rotary balance, it should be pointed 
out that proper settings of the instantaneous spin radius (distance between 
the Instantaneous axis of the rotation and the center of gravity of the 
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model) arid the headiiig or azimuth angle must he Included in order to 
simulate fully (exc^t for acceleration effects) the motion of a free 
airplane model in a steady spin or at any instant in the recovery, ^y 
using the correct values of the instantaneous spin radius , the azimuth 
angle ^ the angles hetvreen the model hody axes and the instantaneous spin 
axis, fl-ud. the tunnel airspeed and hy rotating the model at the proper 
nn gni rt velocity, all the variables of the motion (that is, u, v, w, 
p, q., and r) may he obtained simultaneously. At the time the aero- 
dynamic data used in this Investigation were accumulated, it appeared as 
though an impracticably large amount of testing time would be required 
if the spin radius and heading angle of the model were included as vari- 
ables in obtaining the data for interpolation; however, previous Tinpub- 
llshed work shows that, for apparent steady-spin conditions, rotary- 
balance measurements with zero spin radius did not differ appreciably 
from measTirements where spin radius was used when a, p, and ft were 
adjusted to be the same for the two conditions. If the radius of turn 
should become relatively very large, however, as it may near the end of a 
recovery, the use of zero-spln-reidius rotary-balance data in calculating 
the spin recovery may have some effect on that part of the calculated 
motion. 

In a previous investigation (ref. 7) which the rotary balance 
was Tised to measure the aerodynamic forces and moments under conditions 
which were representative of apparently steady spins in the free-spinning 
tunnel, even though spin radius was used, the values obtained for the 
aerodynamic forces and moments were fotmd to be Inconsistent with the 
assimption of a steady spin in that the aerodynamic and inertial forces 
and moments did not balance. Because of this discrepancy either the 
aerodynamic data or the inertial characteristics must be modified to 
simulate equilibriian in the steady spinning condition. For the present 
investigation also, the inertial and aerodynamic forces and moments did 
not balance; thus, the aerodynamic data were arbitrarily modified by 
adding the proper Increment to each force and moment coefficient to obtain 
equilibrium in the steady spin. During the calculated recovery motion, 
the aerodynamic data were modified by adding these same increments at 
each Instant. The use of these corrections does not mean that the balance 
gives erroneous readings, and as mentioned previously, the inconsistencies 
are believed to be due to difficulties involved in mounting the model. 

The Increments added to the force and moment coefficients to obtain 
equilibrl-um in the free-spinning condition at t = 0 were evaluated as 
follows: !Ehe rudder-neutial aerodynamic coefficients for the free- 

spinning condition were obtained from the tabulated nidder-neutral rotary- 
balance measurements by Interpolating for the proper values of a, p, 
and fl. Since the free-spinning condition had occurred when the rudder 
was deflected in the pro-spin direction, the proper rudder-with-spin 
Increment was then added to each aerodynamic coefficient. This increment 
was obtained from the tabvilated rotary-balance measurements of mdder-with- 
spin increments by interpolating for the proper values of a and fl. The 
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rudder-deflection increment^ were assumed to “be constant with The 

necessary equilibrium aerodynamic increments were obtained hy subtracting 
the rudder-with-spin a^odynamic coefficients from the free-spinning 
model inertial coefficients. 

At t = 0 the spin-recovery motion was initiated hy throwing the 
rudder against the spin. In the equations of motion, this initiation of 
the spin recovery was accon^lished hy subtracting the rudder-with-spin 
increments from the equilibrium aerodynamic coefficients and then adding 
the rudder-eigainst-spin increments. The values for the attitudes, veloc- 
ities, and rudder-against-spln aerodynamic-coefficient factors when time 
equsQjed zero were substituted into the equations of motion along with the 
mass factors of table I, and initial' values of u, v, ,w, p, 4# ^ 

were calculated. These values were assumed to remain constant over a 
small increment of time, and new values of the velocity terms at the end 
of the small time increment were calculated by standard acceleration- 
velocity relationships such as 


^ 

and 

Pn = Pn-1 + 


New valvies of the attitude terms in the equations of motion at the end 
of the small time increment were calculated by the relationships 


(- sin 9e)^ =' (- sin + 


d(- sin 0g) 
dt 


At 


(cos 9e sin = (cos 9g sin 0e)n-l 


d(cos 9g sin 0g) 
dt 


At 


and 


cos = (cos 9^ 


Tdiere the last term in each relationship was assumed to re main constant 
over the small time increment. Values for these last terms were 
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determined "by the method presented in appendix A. A new value of Vj^ was 
obtained hy calculating the vectorial sum of the new values of u, v, and 
w. New valties of the aerodynamic force and moment coefficients at the end 
of the small time interval were obtained from the rotary-balance data in a 
manner similar to that discussed previously for obtaining the coefficients 
at zero time; namely, the same infcerpolative procedures were used in the 
rudder-neutral data, increments were added to the coefficients to account 
for the rudder being against the spin, and the same correction factors 
(arbitrarily assumed to be constant correction factors to the balance 
data) alrea^ discussed as being necessary for equilibrium were added. 

37 he new values of a and 3 needed for Tnalcing the interpolations were 
determined by ttsing the new values of linear velocities in the formulas 
noted previously which relate u, v, w, and Vg to a and p. !The 
new value of fl needed for making the interpolations was taken as the 
vectorial stna of the new valties of p, q, and r. 

The step-calculation process was repeated many times to obtain 
time histories of the individual attitude, velocity, and acceleration 
factors In the equations of motion by svibstltuting the newly obtained 
values at each step In the right-hand sides of the equations of motion. 
Time-history plots were made of these variables and also of a, p, 
and fl. A time-history curve Indicating the number of spinning turns 
the model had made in space ^ at any time was obtained from a graphical 
integration of a curve of ijr against time which was calculated by 

if 

e sin 0Q 

a relationship mentioned prevlovisly. In the grapMcal-integration tech- 
nique, the area under a portion of the curve of ijfg against time for a 
given small time increment represented the change occurring in the value 
of ilTg dioring that time Interval. 

Time histories for the angles ^ and 0 ’ were calculated from 

0 = sln”^ (sin jZig cos 0e) (fig* 3) 


and 




sin 


-1 




The attitudes of the resultant force vector acting on the airplane, with 
respect to the XZ and XT body planes, were calculated by (see fig. Ij-) 
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= sin" 


V + ur - wp 


r(u + wq. - vr)^ + (v + tjt - yp)^ 


/F2 = sin” 


w + vp - uq. 


/(u + wq. - vr)^ + (v + nr - wq)^ + (w + vp - uq)^ 


and the attitudes of tbe axis of resultant rotation vith respect to these 
planes were calculated "by (see fig. 5) 


= sin 


2 2 

? + q. 


Tjj = sin-^ ^ 


The angle a was determined l3y 


cos a s= cos Tq cos Ty^ cos ^6jj - + sin sin Ty^ 

(fig. 6) and the velocity of the airplane along the axis of resultant 
rotation was calculated hy 


V = Vj^ cos a 


The inclination of the axis of resultant rotation from the vertical at 
various times d\n:ing the recovery was determined hy (appendix B) 


cos Kjj = cos tti cos + sin sin cos IH 


The velocity of descent toward the earth at various times during the 
recovery wets obtained hy adding vertical velocity components, as 
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Vy = - 11 sin 0 q + v cos 0^ sin + w cos 0g cos 

and the inclination of the resultant velocity vector from the vertical 
■was calculated by 


cos Krr 


VV 

Ve 


The calculations -were continued until the plotted time-history 
curves indicated that the airplane had achieved such an attit-ude and 
motion that a recovery from the erect spin had occurred. 


RESULTS AMD DISCUSSION 


The calculated motion is presented as time-history plots of the 
varioTis attitijde, velocity, and acceleration factors in figures 7 "to 15* 

Becaiise some of the modifications to the aerodynamic data -were rela- 
tively large, the detailed calculated recovery motions nay not r present 
those of the model actually tested but may be more representative of the 
recovery motion of a hypothetical airplane. The adjustments made to the 
aerodynamic data, however, did not modify the slopes of the aerodynamic 
data with respect to the variables of the motion and the modifications 
might be considered as changing the aerodynamic characteristics of the 
model in order to obtain tidm in the observed steady spinj hence, the 
stability of the airplane in the spin was not affected. The increments 
in the meas\n:ed aerodynamic coefficients arising from the rudder deflec- 
tion and from the changes in the variables of the motion during the 
recovery do, of course, change the slopes of these aerodynamic coeffi- 
cients. It is believed that these latter increments and the inertia 
characteristics of the model were obtained -with good accuracy. There- 
fore, in the event that the stability of the airplane in the spin deter- 
mines the more important characteristics of the recovery motion, the 
calculated recovery may fairly well r^resent that of the actual rotary- 
balance model tested. 

As may be seen from the figtires, the plots show initial gentle oscilla- 
tions from the original, values of the attitude, velocity, and acceleration 
factors, -with the oscillations gradually Increasing in magnitude -vd-th elapsed 
time. U^ to a time of 5 seconds, the changes that occurred in the various 
factors -were relatively small. During the sixth second, large changes 
occurred. At the end of 6 seconds, had reached a -value of -97° (fig. 7), 

the angle of attack had reached zero ffom its original hi gh positive value 
and -was going negative (fig. 8), and -ijig had almost gone to zero (fig. 9). 
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These factors and their rates of change indicate that at 6 seconds the air- 
plane had made a rather abrupt left-wing-down rolling motion to a quasi- 
inverted attitude in space and was still rolling, and recovery from the 
original spin appeared completed. A pilot in the corresponding aii^lane 
should at about this time readjust his controls, probably to neutral or 
near neutral, to achieve a dive after the recovery motion and to prevent 
the airplane ftom entering another developed spin either inverted or 
erect. The delayed-action, rather abrupt recovery motion indicated by 
the calculations appears to have its counterpart in the recoveries from 
spins of some full-scale airplanes. For these airplane recoveries, the 
controls applied for recovery must be maintained by the pilot for an 
appreciable amount of time before the spinning motion is affected greatly, 
after which the recovery occvlps rapidly. The rather abrtqat left-wing- 
down rolling motion indicated by the calculations occurred just after 

(fig* 7) bad reached its hipest positive value (rigjht wing down) 
during a phase of the airplane oscillations. From analysis of the 
rotary-balance data, it appears that this rolling motion was the result 
of a large negative rolling moment (probably due to the phenomenon of 
effective dihedral) which built on the airplane as the inward (posi- 
tive) sideslip reached a maximum (fig. 8). 

An interesting feature of the cxirve (fig. 9) is the sudden 
large increase which occurred at 5*8 seconds, just before ijre decreased 
rapidly to almost zero. This increase which represents a speeding up of 
the rotation about a vertical axis was probably associated with the phe- 
nomenon of conservation of angular momentum in that it occurred as the 
airplane pitched down to a steep attitude in space (Se^ fig. 8) and thus 
decreased the moment of inertia about the vertical axis. 

As may be seen from figure 11, the resultant linear velocity Vr 
increased subtly above its original value during the recovery. Tbrou^- 
out the recovery motion, the veloclly of descent toward the earth Vy 
(fig. 15 ) remained about the same as the resultant velocity Vr, an indi- 
cation that the path of the airplane center of gravity and the direction 
of the resultant wind remained almost vertical. This indication that the 

resultant wind remained almost vertical is also reflected in the kw 

''R 

curve in figure 15 . Figure 12 shows that the resultant force was directed 
above the XT body plane, forward of the TZ body plane, and alternately to 
the left and rigpit of the XZ body plane throu^out the calculated motion. 
Figure 15 indicates that during the spin and most of the recovery, the 
axis of resultant rotation remained almost parallel to the plane of symme- 
tiy of the airplane as evidenced by the fact that bQ remained within a 
few degrees of zero for the first 5 seconds after rudder reversal. Only 
during the last second of recovery did the axis of resultant rotation 
deviate appreciably from the plane of symmetry. The angle a between 
the resTiltant linear velocity and the axis of resultant rotation increased 
a great amount as the motion progressed towEtrd a recovery (fig. l^l-). 
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Along with this increase in cr, the velocity component of the airplane 
parallel to the axis of result^t rotation decreased correspondingly 
(fig. 15). As V decreased, the angle between the axis of resultant 

rotation and the vertical increased until, at the time of the roll-over 
recovery, the axis of resultant rotation was almost horizontal in space 

(fig. 15 ). 

The attitude and motion of the airplane at any time during the 
recovery may be rather completely discerned by slTmltaneous reference to 
the time-history conves. 

As may be seen from figure 9 , the calculations indicated that the 
airplane made about 2 turns in space while in the recovery motion. Of 
interest is the fact that the results obtained with the free -spinning- 
tunnel model mentioned previously also showed an approximately 2-tum 
recovery. This agreement indicates that the calculation method may give 
qualitative accuracy in spite of the presently unexplainable equilibrium 
discrepancies. Any application of this method to obtain the details of 
a spin recovery for a specific airplane depends on the possibilily of 
explaining these equilibrium discrepancies. 


CONCLUDING EMAEKS 


A step-by-step calctilation has been made of an ai3T>lane spin-recovery 
motion by use of modified wind-tunnel rotary-balance measurements, 
applicable equations of motion, and spin-geometry relationships. Diffi- 
culties encountered in applying the rotary-balance data in the calcula- 
tions are discussed, and it is pointed out that certain inconsistencies 
must be cleared up before the method can be accepted as adequate to give 
detailed spin-recovery motions for a specific airplane. Time-history 
plots are presented which represent variations, after rudder reversal for 
attempted recovery, of ensiling attitudes, velocities, and accelerations of 
the airplane with resjiect to the airplane body axes, the earth, the rela- 
tive wind, and the resultant axis of the spinning motion. The time- 
history curves indicated initiaJ. small oscillations of the factors from 
their original values, with these oscillations increasing gradually in 
amplitude and with nearly all the lerge significant changes occurring 
near the end of the time required for recovery. 


La n gley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., March 5, 195^. 
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APPENDIK A 

METHOD OF DETEEMUSTENG ATTITUDE TEEMS 


As noted in the body of the paper, escpressions for angular velocities 
about principal axes in terms of Euler’s space angles are 

P = 9^e “ ife ®e 
q. = 0e cos + iifg cos 0^ sin 
r = -0Q sin cos 0g cos 

Eatpressing these relationships in terms of as 

. $le - p q. - Qe cos 0g r + 0g sin 0g 
sin 0g ” cos 0g sin 0g ~ cos 0g cos J^g 

and cross-multiplying the resulting equalities, two at a time, gives 

cos 0g sin 9^e(?^e “ P) “ 0e(*l ~ ®e ^e^ = 0 

cos 0g cos 0g(q. - 9g cos j5g) - cos 0g sin 5^e(r + 0g sin J^g) = 0 

and 

cos 0g cos ^fgCjie - P) - sin 0g(r + 0g sin ^g) = 0 

When escpanded and rearranged these equations become 

• ^ 

(cos 0g sin + sin 0g cos 0g9g) - p cos 0g sin 0g - q sin 0g = 0 

(- cos 9e9g) + q cos 0g cos ^g - r cos 0g sin ^g = 0 



NACA TK 3188 


19 


and 

(cob 0e cos 9fgJ!(g - sin 0Q sin - P cos 0g cos - r sin 0 q = 0 

Replacement of the terms in parentheses hy their respective equivalents 

- 0e cos 

A(. 3in 9e) 

and 

^(cos 0e sin 0g) 

gives 

~(cos 0e cos ^e) = ~ P cos 0e sin " 9. sin 0e 

CL*C 

— (- sin 0e) = - q. cos 0e cos + t cob 0g sin 


and 


—(cos 0e sin 0e) = P cos 0© cos + r sin 0e 
dt 


The latter three equations were xised to determine values of the left- 
hand sides at zero time and also at the beginning of each time interveil 
hy STibstituting ang nift-r velociiy and attitude values at correspon ding 
times in the ri^t-hand sides of the equations. 
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APPENDIX B 

MEEHOD OP DECEEMINING INCLINAfPiaN OF AXIS OF RESULTANT 
ROTATION FROM VERTICAL 


The inclination of the axis of resultant rotation from the vertical 
was determined hy using spherical trigonometry and by making use of the 
values calcvQated herein for the Inclinations of the axis of resultant 
rotation with respect to the airplane and for the attitudes of the air- 
plane with respect to the vertical. The method is illustrated by the 
following sketch and equations -vdilch indicate the spherical trigonomet- 
rical procedure used for a condition in which the airplane had positive 
values of p, q., and r and was in an erect attitude with its right 
wing somewhat below the horizontal (^e is positive, < 90 °) • Imagine a 
sphere of arbitrary radius having its center at the airplane center of 
gravity. Let points labeled X, Z, ft, and Vy in the sketch repre- 
sent, respectively, points where the surface of the sphere is inter- 
sected by the positive directions of the X and Z body axes, the axis of 
resultant rotation, and a verticeil line through the center of gravity. 
The spherical triangles indicated are then created. 


X 



As may be seen from the appropriate spherical triangles, the m ag ni tude 
of angle I may be obtained from 


cos (90^ - 05) = cos 90° cos + sin 90 ^^ sin ct^ cos I 
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in which 


sin 


cos I = 


“5 


sin 

the magnitude of angle II may he obtained from 


cos 


/ -1 
(cos — 

\ “ 


= cos /cos ^ ^cos 90° + sin 90 ° sin (sin ^ Vr^ + q £ 

\ «/ V « > 


in which 


cos II 


E = 

a a 


cos n = 


Vr2 + q2 


and the magnitude of angle III from 


in = H - I 


Ihen^ an the appropriate spherical triangle in the sketch shows ^ 
cos Kjj = cos a£ cos ogL + sin a£ sin cos III 


In ■using this method care -was exercised so that a value of angle III of 
proper magnitvide -vreis obtained for each instant of time at which the 
angle kq 'VTas calculated. When the vectors representing a vertical line 
through the center of gravity and the axis of resultant rotation were on 
opposite sides of the XZ body plane, a plus sign -was required in the for- 
mula relating angles IH, H, and I to obtain the proper magnitude of 
angle TTT . For other possible locations of these -vectors not encountered 
in the present calcvilated motion, supplementary angles of I and H or 
both may have to be considered. 
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TABLE I — MASS CHARACTERISTICS, CONTROL SETTINGS, AND 
SPIN CHARACTERISTICS FOR AIRPLANE CONFIGURATION 


Mass characteristics: 
Wei^t, lb 

c 

c 

p. at 15,000 ft altit\:ide 

ix 

lY 

iz 




mb^ 


^Y - IZ 

mb^ 

^Z ~ % 

wib2 


• • 17,835 
. . 0.212 

. . 0 .009 
. . 17.35 

. . lT,3i^2 

• . 37,920 
. . 53,396 

-14T X 10"^ 

-no X 10“^ 
257 X 


Control settings, deg: 

Elevator, up (stick back) 20 

Ailerons, against spin (stick left in spin to pilot's right) . . . l4 

Rvidder with spin (right pedal forward in spin to pilot's 

right) 50 

Rvidder reversed for recovery to against spin 30 

Spin characteristics: 

fg, radlans/sec 2.165 

9g, radlans/sec 0 

^g, raxllans/sec 0 

0 g, deg -4k 

O, deg 0.4 

Vy, ft/sec 216 

R, ft 6.62 

a, deg 46 

P, <3fig -3.^ 
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(a) 0g and = 0. 




(c) 0g and ilTg =s 0, and in -this case 0 = 0g, 
Figure 1.- Body system of axes and related angles. 
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XZ body plane ^ 


Vertical plane through Y body axis 
IIHIIIII Vertical plane through Z body axis 
1 — I YZ body plane 
I I Vertical plane through X body axis 

Figure 3 »- Sketch indicating spherical trigonometry relationships pertinent 

to determination of angles and o«. Sin ^ . 

^ cos 0 Q 
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Figure 5.- Vector diagram indicating relationships for obtaining the 

angles 6 q and Tq . 
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Time , sec 


Figure 7*“ Calculated variations of 0, and 0* with time during 

recovery motion. 
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Figure 8 — Calcsulated variations of Og, a, and ^ with time during 

recovery motion. 
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Figure 9 »- Calculated variations of turns of model conpleted in 

space with time during recovery motion. 
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Figure 10.- Calculated variations of angular velocities and accelerations 

■with time during recovery motion. 
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Figjure 11,- Calculated ■variations of linear -veloci-ties and accelerations 
, . -with time during recoveiy motion. 
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Time , sec 


Figure 12.- Calculated variations of Z.Fi and 1^2 with time during 

recovery motion. 
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Figure l4.- Calcvilated variation of a with time during recovery motion. 


ft/sec 


38 


KACA TN 5188 




Figure 15.- 


Calculated 


variations of Vy, V, 
during recovery motion. 


and 


Kjj with time 


NACA-Langle; - 6-18-84 - 1000 


